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Abstract — Recently, a method for obtaining diffraction-attenuation resistant beams in 

absorbing media was developed through suitable superposition of ideal zero-order Bessel 

beams. In this work, we will show that such beams maintain their resistance to diffraction 

and absorption even when generated by finite apertures. Also, we shall extend the orig 
O 

inal method to allow a higher control over the transverse intensity profile of the beams. 

Although the method has been developed for scalar fields, it can be applied to paraxial 



vector wave fields as well. These new beams can possess potential applications, such as 
free space optics, medical apparatuses, remote sensing, optical tweezers, etc.. 
OCIS codes: 050.1940; 070.2580; 090.1970; 140.3300; 170.4520; 230.0230; 260.1960; 
350.7420; 999.9999. 

1 Introduction 

About six years agofTJ, |2], an interesting method was developed, capable of delivering 
beams (in non-absorbing media) whose longitudinal intensity pattern (LIP) could be 
previously chosen. This method, named "Frozen Waves", is based on the superposition 
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of co-propagating Bessel beams, all of them with the same frequency. 

A few time later, this method was further generalized [3], allowing us to model the LIP 
of propagating beams in absorbing media. As a particular case, diffraction-attenuation 
resistant beams were obtained, that is, beams capable of maintaining both, the size and 
the intensity of their central spots for long distances compared to ordinary beams. 

The referred method for absorbing media was developed from appropriate superposi- 
tions of ideal zero-order Bessel beams. This has two fundamental implications: a) beams 
with an infinite power flux (due to the use of ideal Bessel beams), and b) beams with a 
spot-like transversal profile (due to the use of zero-order Bessel beams). 

In this paper, we shall extend the above method in such a way as to obtain a more effi- 
cient control over the transversal intensity profile of such beams by adopting superposition 
of higher order Bessel beams. 

We will also show that diffraction-attenuation resistant beams can maintain their 
interesting characteristics even when generated by finite apertures, i.e., even when we 
transversally truncate the Bessel beams that composes the desired beams. This demon- 
strates that, once the generation scheme is chosen, be it with antennas (in microwaves 
and millimetric waves), holograms or spatial light modulators (in optics), the resultant 
beam may possess characteristics analogous to those of the ideal case, at least until a 
certain field depth. 

Although the method has been developed for scalar fields, it can be applied to elec- 
tromagnetic waves in the paraxial regime and we will elucidate this point. 

These new beam solutions can possess potential applications in medicine, remote 
sensing, free space optics, optical tweezers, etc.. 

2 Resume of the method: Ideal Diffraction- 
Attenuation resistant beams in absorbing media 

In [3] one of the authors developed a method capable of furnishing, in absorbing media, 
beams that are resistant to the effects of the diffraction and, the most important, capable 
of assuming any previously chosen LIP, in p = 0, in the range < z < L. 



As a particular case, this method is capable of furnishing diffraction-attenuation resis- 
tant beams for long distances compared to ordinary beams. With " diffraction-attenuation 
resistant" we mean that the central spots of these new beams maintain their shapes and 
also their intensities while propagating along an absorbing medium. 

Such a method is based on suitable superpositions of equal- frequency zero-order Bessel 
beams and the core idea is to take on extreme the self-reconstructing property of the 
non-diffracting waves [I]- [28]. The obtained beams possess (initial) transversal field dis- 
tributions capable of reconstruct not only the shapes of the central spots, but also the 
intensities of these spots. This happens without active action of the medium, which con- 
tinues absorbing energy in the same way. This section presents the method developed 
in [3] with more details and some news. 

In an absorbing media with a complex refractive index given by 

n{uj) = tir(u)) + iriiiuj) (1) 

a zero-order Bessel beam can be written as 

*P = Jo [{k P R + k pI )p] exp(i/3 R z)exp(-/3 I z)exp(-iut) (2) 

where /3r + ifii = j3 and k p n + k p i = k p are the (complex) longitudinal and transversal 
wave numbers, respectively, with their real and imaginary parts being given by: 



Pr = cos u ; pi — cos 6 (S) 

c c 

and 



n R u . mu . ( 

k pR = sin 9 ; k pI = sin 9 , (4) 

being < 9 < ti/2 the axicon angle of the beam. Notice that k 2 = n 2 u 2 /c 2 — 1 . 

One can clearly see that the Bessel beam (T5]) suffers an exponential decay along the 
propagation direction "z" , due to the term exp(— fiiz). 



The absorption coefficient of a Bessel beam with an axicon angle 9 is given by 
a.Q = 2/3/ = 2n I ucos9/c, its penetration depth being 5g = 1/ag = c/(2corijcos9). It 
is interesting to notice that, because the transverse wave number k p is complex, the beam 
transverse profile decays as a Bessel function until p ~ l/2k p j, beyond which it will suffer 
an exponential growth. This physically undesirable behavior occurs because Eq.(j2]) rep- 
resents an ideal Bessel beam which needs to be generated by an infinite aperture. This 
problem, however, is solved when the beam is transversally truncated, i.e., when we use 
finite apertures to its generation. In these cases, the exponential growth along the trans- 
verse direction (for p > l/2k p i) must cease for a given value of p, and when the radius R 
of this aperture is such that R < l/2k p j, this exponential growth does not even occurs[3]. 

It is important to remember that the efficient generation of a Bessel beam[5j[6] occurs 
when the aperture radiusj is such that R >> 2A/k P R. 

From the two conditions for R mentioned above, one can showj3] that, in an absorbing 
medium, a Bessel beam generated by a finite aperture of radius R will possess acceptable 
characteristics when ur » nj, i.e., when the coefficient of absorption is such that a << 
l/\->6»\. 

All cases considered here must obey this condition. 

Now that the basic characteristics of a Bessel beam in an absorbing medium are 
understood, let us present the method developed in[3]. 

The idea is to achieve, in an absorbing medium with a refractive index n(u) = 
n i?( w ) +inj(u) , an axially symmetric beano, \l/(p, z, t), whose LIP along the propagating 
axis (i.e., on p = 0) can be freely chosen in a range < z < L. Let us say that the desired 
intensity profile in this range is given by |F(^)| 2 . In order to obtain a beam with such 
characteristics, the following solution is proposed: 



*When generated by a finite aperture of radius R>> 2A/k P R situated on the plane z = 0, the solution 
@ becomes a valid approximation only in the spatial region < z < R/ tan 9 = Z and to p < (1 — z/Z)R 
tin this paper we use cylindrical coordinates (p, </>, z). 



with 



N 

9(p,z,t) = £ A m J (k PmP )e i ^e~^ t 

m=-N 

(5) 

N 

= e~ iut J2 A m J ((k pRm +ik pIm )p) e^^e-^ , 

m=—N 



& ="*£-&, (6) 



being 



0Rm _ _ k pR m _ . ^R 



(7) 



h plr, 

where (3 m = /3 Rm + if3 Im and k Pm = k pRm + ik pIm . 

Equation (JSJ) is a superposition of 2 A^ + 1 co-propagating Bessel beams with the same 
angular frequency uj. In (j5]), the coefficients A m , the longitudinal (/3 m ) and the transverse 
ik pm ) wave numbers are yet to be determined. The choice of these values is made such 
that the desired result (i.e., |\]/(p = 0, z,t)\ 2 = |F(z)| 2 within < z < L) is obtained. 

The following choice is made: 



_, 2irm ,„ s 

P Rm = Q + -r- , (8) 



where Q is a constant such that 



^ 27rm a; , s 

< Q + — r < n R - , (9) 

for -JV < m < N. 

Condition ([9]) ensures forward propagation only, with no evanescent waves. In (JH]) 
Q is a constant value that can be freely chosen, as long as (Q be obeyed and it plays 
an important role in determining the spot size of the resulting beam (lower values of Q 
implies in narrower spots) as we will show soon. Besides, Q can be chosen such as to 
guarantee the paraxial regime in an electromagnetic beam. In this case, Q must possess 
a value close to n R u/c. We will see this in details in section 3. 



With the choice (jSj), the superposition fl5]) is written as: 



N 

■ 27T771 . 



*(p,z,t) = e- i » t JQ* £ A m J ((k pRm + ik pIm )p) e^ 2 e- ft ™ 2 , (10) 

m=-N 

where, by inserting (jSj) into (J7J), 

_ /_ 27rm\ m 

V L J n R 

and /c Pm = /c pflm + ifc p / m is obtained through Eq.©. 

The maxima and minima of the imaginary parts of the various /3j m are given by 
(/3/)mm = (Q — 27iN/L)rii/riR and (/3/) max = (Q + 27rN/L)ni/riji, and the central value 
(for m = 0) being given by (/3/) m=0 = Qn I /n R = f3 T . 

Now, let us consider the ratio 



j 



^ _ (/3j)max - (/3j)min _ , iV" /^x 

For A << 1, there are no considerable numerical differences among the various /3j m 
and we can safely approximate them by (3i m w /3 7 , which implies that exp(— /3/ m z) « 
exp(— /3 7 2;). In these cases, the series in Eq.flH]) evaluated on p = can be (approximately) 
considered a truncated Fourier series, multiplied by the function exp(— fijz). Therefore, 
this series can be used to reproduce the desired longitudinal profile |F(2;)| 2 (on p = 0), 
within < z < L, as long as we make 

A m = \ [ L F(z)e^ z e- i2 -^ z dz (13) 

L Jo 

Essentially, this is the method developed in [3J. 

It is interesting to note that countless beams with the same desired LIP, but with 
different values of the parameter Q can be constructed. The basic difference among them 
will be their spot radius (Ap), which can be estimated as being 



2.4 2.4 , , 

Ap *-k = / 2 2 (14) 



So, besides choosing the desired LIP of the beam in an absorbing medium, we can 
stipulate its spot size. In section 4 we will show that a more efficient control over the 
beam transverse intensity pattern is obtained by using higher order Bessel beams in 
superposition (jSJ). 

In resume, we wish to obtain a propagating beam in an absorbing medium possessing, 
inside the interval < z < L, a previously chosen LIP (on p = 0) given by |F(z)| 2 . 
To achieve such a profile, we write the desired beam as a superposition of 2N + 1 co- 
propagating Bessel beams, Eq.flS]), by suitably choosing the longitudinal (j3 m = /Si? m +i/3j m ) 
and transverse {k Pm = k P R m + ik p j m ) wave numbers, and the coefficients A m according 
to Eqs. (l8f9jlllf6jll3p . We can also stipulate the spot radius, Ap, of the beam through a 
suitable choice of the parameter Q in Eq. (Tl4"]) . 

The method demonstrates to be efficient in situations whereja << 2/ A and A = 
47rN/LQ << 1, allowing us to obtain a great variety of beams with potentially interesting 
intensity profiles as, for example, beams capable of maintaining not only the size of their 
central spots, but also the intensity of these spots until a certain chosen distance along 
the absorbing medium in question. We can call these types of beams as "diffraction- 
attenuation resistant beams". 

We finish this section with an exampldll of the above mentioned beam. 
Example 

Consider an absorbing medium with a refractive index n = 1.5 + i0.49 x 10 -6 in 
A = 308nm (i.e., u = 6.12 x 10 15 Hz). In this bulk, at this angular frequency, we have 
the following behaviors for the following wave solutions: a) A plane wave possesses a 
penetration depth of 5 = 1/ot = c/(2umi) = 5 cm ; b) A gaussian beam with initial spot of 
radius 5.6 fim., besides suffering attenuation, is also affected by a strong diffraction, having 
a diffraction length of only 0.6 mm; c)An ideal Bessel beam with central spot of radius 
5.6 /xm (which implies in an axicon angle of 0.0141rad) can maintain its spot size, but 
suffer attenuation, possessing a penetration depth of 5g = 1/ag = c/^urij cos 9) = 5 cm. 



■'■Fortunately, these conditions are satisfied for a great number of situations 
§The same that was given in [3]. 



Now, we are going to use the method exposed in this section to obtain a beam with 
spot radius Ap = 5.6 /tin, capable of maintaining the size and the intensity of its central 
spot until a distance of 25 cm, i.e., a penetration depth 5 times greater than those of the 
Bessel beam and of the plane wavqj, and a diffractionless distance hundred times greater 
than that of the gaussian beam. We also demand that the spot intensity suffers a strong 
fall after the distance z = 25 cm. 

This diffraction attenuation resistant beam can be obtained through solution (|3J) by 
choosing the desired LIP |F(;z)| 2 (on p = 0), within < z < L, according to 

F(z) = { l f ° r ° " Z ~ Z (15) 

(0 elsewhere, 

putting Z = 25 cm, with, for example, L = 33 cm. The stipulated spot radius, Ap = 

5.6 /im, is obtained by putting Q = 0.9999n^w/c in Eg. ( CHI) . This value of Q is used 

for the 0r in Eq.®, ancUJI we choose iV = 20. Finally we use Eg. ({TBI to find out the 

coefficientO A m of ( ITUl) . defining in this way the desired beam. 

We can see that the resulting beam fits very well all the desired characteristics, as is 
shown in Fig.l, being 1(a) the 3D field-intensity and 1(b) the orthogonal projection in 
logaritmic scale. 

The resulting beam possesses a spot radius of 5.6 pm and maintains the size and 
intensity of its central core till the desired distance of 25 cm, suffering after that an 
abrupt intensity fall. 

As we will see in Section 5, the truncated version of this beam will maintain these 
characteristics if the radius R of the finite aperture used for truncation obeys R > 3.8 mm. 
This is already suggested by the Figure l.b of the ideal (i.e., not truncated) beam. 



^In an absorbing medium like this, at a distance of 25 cm these beams would have got their initial 
field-intensity attenuated 148 times. 

"According to ©, the maximum value allowed for N is 158 and we choose to use TV = 20 just for 
simplicity. Of course, using higher values of N we get better results. 

**The analytic calculation of these coefficients is quite simple in this case and their values are not listed 
here, we just use them in Eq. pU]) . 
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Figure 1: (a) 3D field- intensity of the resulting beam, (b) Orthogonal projection of the 
resulting beam in logaritmic scale. 

3 Electromagnetic beams: The paraxial regime 



It should be clear to the reader that the present method is exact, i.e., the obtained 
beams (with the desired LIPs) in absorbing media are exact solutions to the scalar wave 
equation and can possess transverse spots of any sizes, from wavelength dimensions to 
infinity. 

In spite of the method has been developed to scalar fields, it can be used in electro- 
magnetism (optics, microwaves, etc..) in the paraxial regime, where the scalar beam \l/ of 
the previous section would represent the transverse cartesian electric field component of 
a linearly polarized electromagnetic beam, being |\1/| 2 proportional to the time- averaged 
electromagnetic energy density. In these cases, the beam spot size must be much greater 
than the correspondent wavelength. This can be done by choosing the parameter Q of 
Eq.© as Q ps n R u/c. 

We are going to elucidate all these points. 

Consider an absorbing, linear, homogeneous and isotropic medium without boundaries 
and without free charges and free currents. 

The electric (and magnetic) field obeys, in the monochromatic case, the Helmholtz 



equation: 



V 2 E + A; 2 E = 



with E = £(x, y, z)e lu)t and where k is the complex wave number: 



(16) 



k = u^ju(w)e(w) = uJfi(e b + i—j (17) 

being e b (ou) the electric permittivity due to bound electrons, cr{uj) the electric conductivity 
and \x ~ /j,q is the magnetic permeability 
Writing k as 



k = k R + ikj (1< 

and considering the frequency not so close to the resonance regions!], we have [29]: 



k R pau 



ki pa to 



Votb 



yUoCfc 



'1 + 



with the absorption coefficient^ a = 2kj. 
Let E be an electric field given by: 



(-)' + > 



-) - 1 

e b uj 



1/2 



1/2 



(19) 



(20) 



E = E x e x + E z e z , 
and let us apply our scalar method to the cartesian component E x 



N 



E x {p,z,t) = e"^ f ]T A m J (k pmP )e i ^ z 

m=-N 

N 

= e - iu t e i Qz ^ A m J ({k pBm +ik pIm )p)e-^ z e i ^ mz 



m=—N 



(21) 



(22) 



^In this case we can consider both e(,(w) and er(w) real quantities. 

W Not ice that, according to section 2, the absorption coefficient of a Bessel beam is ag = acosf 
2ki cos 9. When 9 — > 0, the Bessel beam tends to a plane wave and ag — > a. 
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where, as before, j3 m = (3 Rm + i(3 Im and k pm = k pRm + ik pIm are chosen according to the 
Eqs. (l8ll9][Tl~1l6"]) . and A m calculated through Eq. (fT3|) . being |F(z)| 2 the desired LIP to the 
electric field component E x , 

Now, to satisfy the Gauss law V ■ D = (no free charges), with E = E x e x + E z e z , 
in the monochromatic regime, we must have: 

E, = -/^* (23) 

Using E x as giving by Eq. (l22|) . we found: 

N k 

E z (p,<f),z,t) = e-*" Yl A m ^J 1 (k pmP ) cos <Pe l ^ z (24) 

m=-N P™ 

The paraxial regime is characterized by beams whose wave vectors of their constituent 
plane waves are almost parallel to their propagation directions. In our case, this direction 
is u +z" and the paraxial regime is reached by putting Q ps tircu/c, implicating that 
\kpml f j3 m \ « 1, which in turn implicates that \E Z \ « \E X \. So, in this circumstance, we 
can use the so called paraxial approximation and write the electric field ( 1211) as 



EnE x e x = (e- iwi ^ A m J (k pm p) e l ^ z I e x , (25) 

\ m=-N J 

(paraxial approximation) 
The associated magnetic field can be found through Faraday's law: 

B = --V x E (26) 

Considering the electric field given by ( l2T|22ll24p , it is easy to show that in the paraxial 
regime the above expression can be approximated as 

B«-i^e y = fe-* wt E A m ^J (k pmP )e^Ae y (27) 

u az v m =- N u J 

(paraxial approximation) 
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Remembering that (3 m = (3 Rm + i(3 Im = (3 Rm + i(n I /n R )(3 Rm , and that < f3 Rm < 
Q + 27im/L, the paraxial regime (Q ~ n R uj/c) implies that (3 Rm ~ n R co/c, and in this 
case /3 m /oj ~ n/c (where n = n r + inj is the complex refractive index). So, we can make 
another approximation in Eq. fl27|) writing: 



B w - ( e-^ 1 Y, A m Mk P mP) e l ^ z ] e y (21 

C » m=~N 



or 



77. ?7 

B k -e, x E = -^e. (29) 

c c 

(paraxial approximation) 
Using Eqs. (j25)29p . we immediately see that the time-averaged energy density for 
monochromatic electromagnetic fields, u = (l/4)Re(E-D* + B-H*), can be approximated 
as 

u « ^Re (e* + i^r J |E X | 2 oc |E X | 2 (30) 

4 \ /j,*c z J 

(paraxial approximation) 
So, we can use our scalar method to obtain, in absorbing media, paraxial electromag- 
netic beams whose the time-averaged energy densities on the propagation axis can assume 
any desired patterns within an interval < z < L. 



4 Extending the method to nonaxially symmetric 
beams: Increasing the control over the transverse 
intensity pattern. 

The method developed in [3] allows a strong control over the LIP (on p = 0) of beams 
propagating in absorbing media. 

Once we have controlled the beam's LIP, the transverse intensity pattern (TIP) can 
be shaped in a limited way; more specifically, the spot size of the resultant beam can be 

12 



chosen by a suitable choice of the parameter Q via Eq.flT] 

In this sectioio we are going to show that it is possible to get a more efficient control 
over the TIP, maintaining, at same time, a strong control over the LIP. It will be possible, 
for instance, to shift the desired LIP from p = to p = p' > 0. In other words, we will be 
able to construct the desired LIP over cylindrical surfaces (instead of over the line p = 0). 
Below we explain this new procedure. 

To obtain these new beams we proceed as before, choosing the desired LIP on p = 
within < z < L, choosing the values of Q and N (observing the Eq.flH])), and calculating 
the values of A m through Eq. ([TBI . 

Having done this, we replace the zero order Bessel beams in superposition ( flOl) with 
higher order ones. The new beam is written as: 



JV 

2tt„ 



tt(p,0,z,t) = e-^VV** £ A m J ll (k pmP )e-^ z e^ mz , (31) 

m=—N 

with p a positive integer and all other parameters (Q, N, L and A m ) given and calculated 
as before. 

For the situations considered here, we have that rij << ur —¥ k p j m << k p R m . This im- 
plies that for each mth Bessel function in Eq. ffTOj) is valid: J^[{k P R m +ik p j m )p] ~ J^(k pRm p) 
for < p << l/k p i m , and in this range each mth Bessel function reaches its maximum 
value at p = p' m , being p' m the first positive root of the equation (d J^{k pRm p)/dp)\ p i m = 0. 

The values of p' m are located around the central value p' m=0 and we can expect a shift 
of the desired LIP from p = to p ~ p' m= o- 

We have confirmed that conjecture in all situations considered by us, mainly in the 
cases where there is no considerable difference among the values of k pRm . 

With this extension of the original method, it is possible to model the LIP over cylin- 
drical surfaces, obtaining very interesting static configuration of the field intensity. In 
particular, we can construct (in absorbing media) hollow beams resistant to the attenua- 
tion and diffraction effects. 



ttThe idea developed in this section generalizes that exposed in Section 5 of reference [2], which was 
addressed to non-absorbing media. 
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To show this method working we are going to obtain a cylindrical light surface of 
constant intensity in an absorbing medium with ur = 1.5 and a = 20m -1 (at A = 308nm 
— y u) = 6.12 x 10 15 Hz), which implies that nj = 0.49 x 10~ 6 . At this angular frequency, 
an ordinary beam would possess a depth of penetration of 5 cm in this medium. 

Let us choose, within < z < L, a LIP as that of section 1, Eq.f fT5|) : 



F(z) = | l ioT °< Z < Z (32 ) 

) elsewhere, 

being Z = 25 cm and L = 33 cm. 

As in section 1, we put Q = 0.9999n R u/c in Eq.flH]). With the chosen values to Q 
and L, the maximum value allowed to N is iV = 154, but for simplicity reasons we choose 
N = 20. 

Using equations ( l8|9lTT|6|l~3|) we evaluate all the (3 m , k pm and A m . But, as we have 
explained in this section, instead of using all these values in ( 1T0|) . we use them in the 
superposition ( 13~TT) . where we choose /z = 4. 

According to the previous discussion, we can expect the desired LIP over a cylindrical 
surface of radius p ~ 5, 318/k P R m=0 = 12.289 /im (that is where the function J-±(k P R m=Q p) 
assumes its maximum value). 

We can see in figure(2.a) the resulting intensity field in a three-dimensional pattern. 
Its orthogonal projection shown in Fig.(2.b) clearly confirms the cylindrical surface of 
light intensity. Figure (2.c) depicts the transverse intensity pattern at z = L/2. It is 
possible to note that the transverse peak intensity is located at p = 12.285 /xm, a value 
near to the predicted one. 

This interesting field configuration is resistant to the attenuation and diffraction ef- 
fects till the distance z = 25cm. We should remember that any other ordinary beam at 
the same frequency, propagating in the same medium, would have a penetration depth of 
only 5 cm. 
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Figure 2: (a) 3D field-intensity of the resulting beam, (b) Its orthogonal projection. 



|T(x,y,z = L / 2)f 




Figure 3: The beam's transverse intensity pattern at z = L/2 

5 Finite aperture generation of the diffraction- 
attenuation resistant beams 



As we have seen, the solution §5§, with 



represents propagating beams in 



absorbing media with the remarkable characteristic of allowing us to choose the desired 
LIP on p = 0, within < z < L, being the spot sizes regulated by the value of the 
parameter Q. The same occurs with solution ( 13~T|) . which in turn allows us to choose the 
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desired LIP on a cylindrical surface. 

Now, we must remember that, in spite of the beams ([5]) and (l3Tj) be exact solutions, 
they do not represent beams generated or truncated by finite apertures. Actually, fields 
given by those solutions in all points of space would require infinite apertures to be 
generated. 

However, if a Bessel beam given by Eq.(J2]) is truncated by a finite aperture of radius 
R » 2A/k p R situated on the plane z = 0, the radiated field - in the spatial region 
< z < R/t&nO = Z (the diffractionless distance of a truncated Bessel beam) and 
< p < (1 — z/Z)R - can be approximately described_!§2] El [6] by Eq.(j2J). 

Taking into account that the solution ( ITU]) is a linear superposition of Bessel beams, 
we can expect that when it is truncated by a finite aperture of radius^] R» 2A/k P R m=N1 
the radiated field - in the regioijJJJ < z < Rj tan# m= _jv = Z m =-iv and < p < 
(1 — z/Z m= _ N )R - will be approximately given by Eq.f fTOj) . 

Now, as we are interested in controlling the LIP of the truncated beam within < 
z < L , we have to guarantee that, after the truncation, all Bessel beams in (TTUj) maintain 
their characteristics until z = L. This is possible if both conditions are satisfied: 



ldd 



^pR m=T 

anc "" 



2 4 
R » t— (33) 



n 2 R u 2 

?2 



(34) 



Rm = -N 



T3 

Z m= _ N > L — > - — > L —$■ R> L, , 

tan6' m= _ A r ^ c 2 pf 

So we can expect that by choosing a finite aperture of radius R large enough to satisfy 
( 133|) and (!34j) . the truncated versions of the diffraction attenuation resistant beams will 
maintain their characteristics, i.e., we will continue to be able to control the beam's LIP. 

In order to confirm our expectation, some examples shall be presented, where we 



^The same is valid for a truncated higher-order Bessel beam. 

^Notice that k P R m=N is the smallest value of all k p n m , therefore if R >> 2A/k p R m=N — > R » 
2A/k p ji m for all to. 

"Here, 9 m is the axicon angle of the mth Bessel beam in ([5]). 
**That is, the shortest diffractionless distance is larger than the distance L. 
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choose some desired LIPs and obtain the correspondent ideal beam solutions, \l/(p, (ft, z, t), 
through Eqs. (!5|8|9lfTTl[6lfT3l) . These ideal solutions, in turn, are used for obtaining 
their truncated versions, ^t{p, 4>, z , t), through numerical calculation of the Rayleigh- 
Sommerfeld diffraction integral for monochromatic waves [30J: 



1 r 2lT r R e ikD ( d \ 

y T ( P ,<P,z,t) = — j o #'y o df/f/—l-^*(j/,<f/,J,t)\ , (35) 

where a circular aperture of radius R, on the plane z' = 0, is used for the truncation, 
being the distance D = J(z — z') + p 2 + p' 2 — 2 pp' cos (0 — 0') the separation between 
the source and observation points. 

A. First Example 

Here we are going to calculate numerically the truncated version of the ideal 
diffraction-attenuation resistant beam obtained in the example of Section II and in ref . [3] . 

The medium in question has refractive index n = n^ + irij = 1.5 + i0.49 x 10~ 6 in 
A = 308nm (i.e., u = 6.12 x 10 15 Hz). The ideal beam in that example was constructed 
to possess (in u = 6.12 x 10 15 Hz) a spot radius Ap = 5.6 pm. and a constant intensity 
of its central spot until a distance of 25 cm. These characteristics were reached through 
the fundamental ideal solution (TTU|) . with Eqs. (181l91[TTf6"][T3"j) . by choosing the desired LIP 
|F(z)| 2 (on p = 0), within < z < L, according to Eq. ffT5|) . putting Z = 25 cm and 
L = 33 cm. The values of Q and iV were chosen to be 0.9999n#u;/c and 20 respectively. 

This ideal beam solution is used as the aperture excitation, ty(p',z',t), in the 
Rayleigh-Sommerfeld diffraction integral, Eq.( l35l) . which is numerically calculated to yield 
\&t(p, z, t), i.e., the truncated version of this beam. 

According to our above discussion, the truncated beam shall possesses a behavior very 
similar to the ideal one provided that the aperture radius R satisfy both conditions (|33| 
13^1) . which furnishes in this case R > 4.9mm. However, due to the fact that the chosen 
on-axis LIP has null value within Z < z < L, we can replace L in condition fl34|) with Z, 
obtaining in this case R > 3.8 mm. We choose R = 3.8 mm. 

After numerical calculation of (|35|) . we obtain the result plotted below. Comparing 
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Figure(3) with Figure(la) of the ideal beam of Section II, we can see that our expecta- 
tions were correct, i.e., by choosing an aperture radius big enough, the truncated version 
becomes very close to the ideal solution in the spatial region of interest. 




p(mm) 



Figure 4: The truncated version of the ideal diffraction-attenuation resistant beam ob- 
tained in the example of Section II 



B. Second Example 

In [3] one of the authors obtained, in an absorbing media, an ideal (i.e., not trun- 
cated) beam presenting an interesting and counterintuitive characteristic. There, it was 
considered a medium with refractive index n = % + irij = 1.5 + z0.46 x 10 -6 (which 
implies in a penetration depth of 5 cm) at u> = 6.12 x 10 15 Hz. At this angular frequency 
an ideal beam was shaped to possess a spot radius Ap = 5.6 /im and a modest exponential 
growth of its intensity until a distance of 25 cm, suffering after this a strong intensity fall. 

In order to reach these characteristics, the desired on axis LIP, |F(z)| 2 , within < 
z < L, was chosen according to: 



F{z) 



(36) 



f exp(z/Z) for < z < Z 
1 elsewhere, 

with Z = 25 cm and L = 33 cm. Taking into account Eq. ffT4|) . it was chosen Q = 
0.9999nf>u/c and the value of N was chosen to be iV = 20. 

With this, the (ideal) desired beam, *&(p, z, t), was obtained in [3] through the fun- 
damental solution (TTUj) . with Eqs. (l8|9lfTTl[6lfT3l) . Below we plotted this resulting ideal (i.e., 
not truncated) beam. 




p (mm) -o.5 



z(m) 



Figure 5: Ideal beam presenting a moderate exponential growth in an absorbing medium 



In this example, we are going to use Eq. (l35j) to obtain the truncated version, 
^t{p, z, t), of the above ideal beam. 

By using conditions ( 1331 13^1) for an efficient finite aperture generation, we obtain that 
the aperture radius should obey R > 4.9 mm. But, for the same reason of the previous 
example, we can adopt R > 3.8 mm, and we choose R = 3.8 mm. 

The numerical calculation of Eq. (|35|) yields the truncated beam plotted below: 

We can see an excellent agreement between the ideal beam and the truncated one 
in the region of interest, confirming that the method works very well in more realistic 
situations close to the experimental ones. 
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Figure 6: Truncated version of the beam with exponential growth 

6 Conclusions 

A few years ago[3] it was developed an interesting theoretical method, from which it 
is possible to construct, in absorbing media, axially symmetric beams whose the LIP can 
be previously chosen. As a particular case, diffraction-attenuation resistant beams were 
obtained, that is, beams capable of maintaining both, the size and the intensity of their 
central spots. 

In this paper we elucidate how to apply this scalar method in paraxial electromagnetic 
waves and extended it to include non-axially symmetric beams, allowing in this way, 
besides the strong control over the longitudinal intensity patter, a certain control over 
the transverse one. With this extension, it is possible to model the LIP over cylindrical 
surfaces. In particular, we can construct (in absorbing media) hollow beams resistant to 
the attenuation and diffraction effects. 

We also use the Rayleigh-Sommerfeld diffraction integral to obtain the truncated 
versions of these new beams, verifying that, provided that the aperture used for truncation 
is big enough, the truncated beams possess the same interesting characteristics of the ideal 
beams. This verification is important because it confirms that the method works very 
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well in more realistic situations close to the experimental ones. 

These new beams can possess potential applications, such as free space optics, medical 
apparatuses, remote sensing, etc.. 
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